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Abstract. In the course of animal development, the shape of tissue emerges in part from mechanical and 
biochemical interactions between cells. Measuring stress in tissue is essential for studying morphogenesis 
and its physical constraints. Experimental measurements of stress reported thus far have been invasive, 
indirect, or local. One theoretical approach is force inference from cell shapes and connectivity, which 
is non-invasive, can provide a space-time map of stress and relies on prefactors. Here, to validate force- 
inference methods, we performed a comparative study of them. Three force-inference methods, which 
differ in their approach of treating indefiniteness in an inverse problem between cell shapes and forces, 
were tested by using two artificial and two experimental data sets. Our results using different datasets 
consistently indicate that our Bayesian force inference, by which cell-junction tensions and cell pressures 
are simultaneously estimated, performs best in terms of accuracy and robustness. Moreover, by measuring 
the stress anisotropy and relaxation, we cross- validated the force inference and the global annular ablation 
of tissue, each of which relies on different prefactors. A practical choice of force-inference methods in 
distinct systems of interest is discussed. 

PACS. 87.17. Rt Cell adhesion and cell mechanics 



1 Introduction 



During tissue morphogenesis, cell- level dynamics, e.g., cell 
morphogenesis, cell rearrangement, cell division, and cell 
death, are orchestrated in time and space to shape the an- 
imal body. As conserved families of signaling pathways in 
morphogenetic processes have been identified, new chal- 
lenges arise, such as how mechanical forces that directly 
act on cells and modify their shapes are integrated with 
biochemical signaling to regulate the correct tissue pat- 
terning 1,23,4.5,6,7. Indeed, a growing number of stud- 
ies now address the mechanical basis of morphogenesis. 
For instance, for epithelial tissue, in which the acto-myosin 
cytoskeleton is connected to a network of cell-cell junctions 
(Fig. HJa)), studies are beginning to clarify how tissue is 
shaped by forces acting along the plane of the adherens 
junction, i.e., tension that shortens a cell contact surface 
and pressure that counteracts the tension to maintain the 
size of a cell (Fig.QJb)) ^^EV^^^M^EMH- 
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Measuring tissue stress is therefore useful for deep- 
ening our understanding of morphogenesis and its physi- 
cal constraints. Various in vivo mechanical measurement 
methods have been developed; these include elastogra- 
phy [18 , photoelasticity |19j . magnetic micromanipulation 
[20] , tonometry [21] , nanoindentation [22] , and monolayer 
stress microscopy (MSM) [53]. Among them, laser abla- 
tion of individual cell junctions is most frequently used as 
a tool to evaluate the tension acting on a contact surface 
of epithelial cells [Tipi]. 

Another approach is based on cell shapes [25l[26ll27l 
|2"5] . If all cells had the same tensions and pressures, all the 
angles between cell contact surfaces would be 120°. Con- 
versely, deviations from 120° would yield information on 
pressures and tensions. If we manage to have the informa- 
tion well posed, we can infer forces and stresses using only 
segmented images (that is, images wherein the cell con- 
tours and vertices have been recognized). Force inference 
is non-invasive; hence, spatio-temporal dynamics of forces 
on more than hundreds of cells can be simultaneously es- 
timated, which represents a distinct advantage over cur- 
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rently available experimental methods. Space-time maps 
of stress obtained by force inference would unveil the novel 
physical principles required to regulate morphogenesis. 

Given the power of force-inference methods, their vali- 
dation in multicellular systems merits thorough and care- 
ful analyses. In the present study, we performed a com- 
prehensive test of force inference. Three types of force- 
inference methods, which differ in their approach to treat- 
ing indefiniteness in the inverse problem between forces 
and cell shapes, were employed (Sect. I2.2.2[) . The first 
method (ST) estimates only tensions, and all the cell pres- 
sures are assumed to be the same. The second method 
(SP) estimates only cell pressures under the assumption 
of uniform tensions. Under these respective assumptions, 
the first two methods treat overdetermined problems with 
respect to unknown variables, that is, cell junction ten- 
sions and cell pressures. The third method (STP) treats 
the ill-conditioned problem and simultaneously estimates 
both tensions and pressures by employing Bayesian statis- 
tics with a prior function representing positive tensions 

EH]. 

Tests of these three force-inference methods were per- 
formed for two artificial data, which was generated by 
numerical simulations. One advantage of using simulated 
data is that we can check the accuracy of estimation by 
directly comparing true and estimated values. One of the 
artificial data is a simulated foam, which constitutes a 
well-studied model system for disordered cellular materi- 
als [33], and the other is a simulated cell population. In 
addition, we also used experimental data from epithelial 
tissues. Patterns of estimated forces were compared among 
different force-inference methods in Drosophila pupal wing 
and notum (Fig.[TJc)). In notum data, we previously intro- 
duced an original type of ablation experiment to measure 
the mechanical state and material properties of a tissue 
[30]. The global ablation method and the force-inference 
methods each depend on different prefactors, and thus the 
comparison between the two methods provides us with an 
opportunity to cross- validate them. Based on the results 
of this study, we will discuss a practical choice of force- 
inference methods for specific purposes. 



2 Methods 

2.1 Systems for tests 

Tests of force inference were performed in two artificial 
systems and two experimental systems. Artificial data are 
generated by simulations of foam and two-dimensional tis- 
sue. The experimental data of Drosophila epithelial tissues 
(wing and notum) used here were reported in our previous 
studies [251150]. 



2.1.1 Numerical simulation of foam 

To create the cluster of bubbles shown in Fig. [2] we use 
the Surface Evolver software [5T] in a mode that describes 



each bubble-bubble interface as an arc of a circle with uni- 
form tension. The Evolver minimizes the following energy 
functional using gradient descent: 

u SE = r kj + J2 Pi ( Ai - A °^ W 

interfaces cells 

where F is the line tension (set to one here), Uj is the 
length of the interface separating bubble i from bubble 
j, and the Lagrange multiplier pi, which ensures that the 
area Ai of each bubble is constrained to an individual tar- 
get value Aq, is the pressure in bubble i. The structure is 
therefore a precise realization of the ideal two-dimensional 
soap froth [2j|3j. 

We start from a polydisperse foam (i.e., a range of cell 
target areas) with 2000 bubbles, a total area equal to one, 
and periodic boundary conditions. This foam is relaxed 
to an energy minimum. The pressure of each bubble and 
the position of each vertex in the circular sample of 244 
bubbles shown are recorded, and since each arc is defined 
by the positions of the vertices at its ends and its mid- 
point, we can calculate the center and radius of curvature 
for each interface. 

In the estimation, bubble-bubble contact surfaces are 
approximated by straight lines, and the validity of the 
approximation is discussed later (Sect. B~51) . 

2.1.2 Numerical simulation of a cell vertex model 

A procedure for generating test data in the cell vertex 
model [12,33 is described in [28 ]. Briefly, the geometry 
of cells is approximated by polygonal tiles with straight 
contact surfaces, each of which is specified by the position 
of vertex = (xi,yi) and its connections. The change 
in cell geometry is determined by relaxing the following 
potential function with Tl processes (reconnection of cell 
contact surfaces) allowed Pl[n][T21pil[T4l[T51[TrJ] : 

= £ f (A.-Avf + £ Aikd + ]T (2) 

* M 1 

Here, = fj — r j indicates the relative positions of the 
ith and jth vertices, and thus |r,j| represents the length 
of the contact surface connecting the ith and jth vertices. 
The first term represents the area elasticity of a cell with 
stiffness modulus K = 100.0 and natural area Aq = 1.2. 
The second term represents the line tension with constant 
value -Tjj, the component of the net tension that is inde- 
pendent of the length of the contact surface. The coeffi- 
cients of line tension, Jy, are randomly selected from a 
Gaussian distribution with mean (r i: ,-) = 0.12 x KAl' 
and variance (AT?-) 1 / 2 = 0.4 x (r i3 ). The third term rep- 
resents cortical elasticity, where Li is the peripheral length 
of the ith cell and the coefficient is set to A = 0.04 x KAq. 

2.1.3 Drosophila wing 

The image collection and analysis of Drosophila pupal 
wings (Fig. [Tfc)) are described in 28]. Briefly, data col- 
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Fig. 1. The structure and force balance of an epithelial 
tissue, (a) Mechanical interactions among epithelial cells act 
mostly in the plane of the adherens junction, where cell ad- 
hesion molecules, cadherin, held cells together. Inside the cell, 
an acto-myosin cable runs along the cell cortex in the plane of 
the adherens junction, (b) Forces acting on a vertex located at 
position ro indicated by a black dot. Tensions along the cell 
contact surfaces pull the vertex in the respective directions in- 
dicated by red arrows, while the cell pressures push the vertex 
in the directions indicated by blue arrows, (c) Two epithelial 
tissues: wing and notum in Drosophila, shown in the adult fly 
(experiments performed in pupa). 



lected in pupal wings at 23 hours after puparium forma- 
tion (h APF) were used for testing the force-inference 
methods in the present analysis. To highlight the shape of 
the cell at the level of the adherens junction, Dacatenin- 
TagRFP was used as a marker of the adherens junction. 
We segmented images by using a custom-made macro and 
plug-ins in ImageJ. We manually corrected the skeletonized 
pattern when necessary. 



2.1.4 Drosophila notum 

For large scale laser ablation, pupae were mounted as de- 
scribed in [3D]. Briefly, we used short laser pulses to sever 
the adherens junctions in an annular region around an 
approximately 30-/im-radius circular tissue domain. The 
wound margin retracted with a speed reflecting the stress 
prevailing in the tissue before ablation, and its anisotropy. 
For each experiment presented in [30], using in- house au- 
tomatic software (M ATLAB [33]) based on a watershed 
algorithm followed by manual checking, we segmented the 
cell contours in the circular patch of tissue before and after 
its periphery was ablated. To increase the image quality, 
the initial image of the tissue at rest was obtained as an 
average of 10 images before ablation. Similarly, the final 
image was obtained as an average of raw images (n = 3 
or n — 1 1 according to the quality of the movie) when the 
ablated circular tissue domain was again at rest. 



2.2 Force-inference methods 

Here, we briefly outline how one infers forces and stress 
from patterns of epithelial cell shapes and their connectiv- 
ity. Detailed descriptions of force inference can be found 
in [2T1I28) . As input, we take a segmented image of ep- 
ithelial cells. Epithelial tissue is approximated as a two- 
dimensional sheet, and cells are represented by polygonal 
tiles. Because the curvature of the cell contact surface is 



small in most epithelia, we approximate a cell contact sur- 
face as a straight edge. Here, the unknowns are the tension 
of each contact surface and the pressure of each cell. If 
the deformation of cells is sufficiently slow in a dissipative 
environment, these forces are almost balanced, and the 
system lies in the vicinity of an equilibrium state (qua- 
sistatic assumption). Then, by considering force-balance 
equations with a given cell geometry, one can infer ten- 
sions and pressures. 



2.2.1 Equations 

This section presents the force-balance equations used to 
deduce forces and stress. Although the derivation of the 
force-balance equations was described in [55], we reiter- 
ated it to make the present paper self-contained. The ge- 
ometry of the tissue is specified by the positions of ver- 
tices, Ti — (xi,yi), and their connectivity. We denote the 
tension of the cell contact surface that connects the ith 
and jth vertices as and the pressure of the ith cell as 
Pi. Let us consider the force acting on the 0th vertex at 
the origin ro = (0, 0) in Fig. QJb). The forces in the x and 
y directions are given by 



»=i 

3 



i=l 



Xi 



(3) 
(4) 



Here, Tj = T^o and P4 = P\. By using the orientation 
of the edge, % = tan -1 (yij/xij), the coefficients of 
are simply expressed as Xij/\rij\ = cos 6^- and yij/\fij\ = 
sin 8ij . Pressures act all along the sides of the cells in their 
normal direction, and projection to the x and y axes gives 
the prefactors — j/j and Xi, respectively, of the normal force. 
Half of the force acts on the end-points of the cell contact 
surface, represented by the second terms in Eqs. ([3]) and 

More rigorously, Eqs. ([3]) and (0]) can be obtained by 
differentiating of a potential function. Consider a potential 
function U({r{\) that determines the tissue's mechanical 
characteristics, as exemplified in Eq. ©. The derivative of 
the potential function with respect to r, gives the forces on 
the ith vertex as Fi = dll({ri})/dri. Fi can be written 
as 



E 



dU d\r ok \ 
d\r jk \ dri 



x ^ dU dAj 

^ OA, dr, 

j J 



(5) 



where \rjk\ and Aj are the length of the cell contact sur- 
faces and the cell area, respectively. With the definition of 
pressure and tension, Pj = —dlf/dAj and Tjk = dll/d\rjk\, 
the above equation leads to Eqs. (J3J and (jl]), irrespective 
of the functional form of U . 

Suppose we have an image in which N cells are sur- 
rounded by R cells. The numbers of cell contact surfaces 
and vertices in the image are denoted as E and V , re- 
spectively. Repeating the same derivation of force-balance 
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equations for every vertex, we obtain a vector F = (F x , F y ) 
that represents the forces acting on vertices in the x and 
y directions as 

F = A T T + A P P = AX. (6) 

Here, T and P are vectors composed of Tij and Pi, re- 
spectively. X = (T, P) represents the unknown variables 
to be inferred. At and Ap (and thus A) are 2V x E and 
2V x (N + R) matrices representing the coefficients of Tij 
and Pi in Eqs. (J3J and (HJ), respectively, and they are deter- 
mined by the positions of the vertices. Under the assump- 
tion of quasistatic cell shape changes, the force-balance 
equation becomes 

A T T + A P P = 0. (7) 

Equation (JJJ gives us a relationship between the observ- 
able geometry (angles and lengths; see Eqs. © and (®) 
of cells and the unobservable tensions T and pressures P. 

In Eq. (0, the scale factor of forces is undetermined, 
because the cell shape does not provide any information 
about it. The force-inference method therefore estimates 
relative values of forces in the scale, as described below. In 
addition, hydrostatic pressure (the baseline value of pres- 
sure) cannot be determined, because Eq. ((7|) is invariant 
under a uniform increase of pressure (see Eq. ([3])). Thus, 
it is differences in pressures among cells that are inferred. 
The estimated tensions and pressures are related to true 
ones as T = cT trU e and P = cP true + Ap. Unless men- 
tioned explicitly, the prefactor c is selected to satisfy the 
requirement that the average of the tensions should be 
unity: 

E W , T «/ S = 1 » ( g ) 

and the hydrostatic pressures are selected such that the 
average of cell pressures is zero: 

X\P<=0. (9) 

As described in our previous study [35J , the critical dif- 
ficulty in calculating forces from Eq. ([7]) originates from 
an insufficient number of force-balance equations to de- 
termine the unique solution of the unknown variables X. 
By considering the topology of a network of cell con- 
tact surfaces, we can show that the number of unknowns 
(N + R + E) is smaller than the number of conditions 
(2V) by R + f + 1, where R and / are the number of 
surrounding cells and the number of four- way junctions, 
respectively. One indefiniteness arises from the hydrostatic 
pressure as described above, and another R + f indefinite- 
ness results from the boundary conditions and four-way 
junctions. The latter indefiniteness should be handled by 
a proper assumption on a system to get plausible esti- 
mates. 

One possible formulation of the inverse problem is to 
decrease the number of unknown variables by assuming a 
relationship among variables (variable-reduction approach). 
For example, Chiou et al. [37] set all cells to have the 



same pressure Pi = Pq. Then, with a constraint of fi- 
nite scale, given by Eq. ((SJ, the inverse problem becomes 
over-determined and has a unique solution. Another for- 
mulation is to adopt Bayesian statistics, which now pro- 
vides a standard framework for this ill-conditioned prob- 
lem, where our expectation for the system is incorporated 
as a prior (Bayesian approach) [35,36]. For example, we 
used a prior function expecting that tension is distributed 
around a positive value |28] . In a Bayesian framework, 
forces are inferred by maximum a posteriori (MAP) es- 
timation after the marginal likelihood is maximized with 
respect to a hyperparameter. 

With the obtained values of tensions and pressure (X = 
(P, T)), one can integrate them to deduce the global stress. 
The stress tensor is evaluated with the Batchelor stress 
tensor given by [25113"?] 

'-if-E^+E^ 1 ^). ex-) 

where I is the two-dimensional identity matrix and A = 
J2i Ai is the total tissue area. Because the scale and the 
hydrostatic pressure are undetermined in the force-balance 
equations, the scale of <x is also undetermined as well as 
the additional pressure —Apl. Some quantities derived 
from the tensor are independent of Ap; they include the 
maximum stress direction and the difference of two eigen- 
values of er. 



2.2.2 Force inference methods to be tested 

In this study, we performed a comparative test on three 
types of force-inference methods, which differ in their ap- 
proach of treating indefiniteness in the inverse problem 
between forces and cell shape. The first two methods re- 
duce the number of unknown variables and treat overde- 
termined problems with respect to unknown variables T or 
P (variable-reduction approach) . The third method treats 
the ill-conditioned problem by employing Bayesian statis- 
tics (Bayesian approach) [28]. We have called them ST, 
SP and STP, where the "S" stands for "straight" edges 
(curvatures are neglected and cells are treated as poly- 
gons) ; "T" and "P" mean that tensions and pressures are 
unknown, respectively. 

Tension inference under an assumption of uniform pressure 
(ST) 

The first method ST estimates only tensions. In ST, the 
difference in pressures among cells is assumed to be small 
and cells are approximated to have the same pressure 
Pi = Pq. Under this assumption, ApPq vanishes, as is 
immediately evident from Eq. ^ and (UJ). Then, Eq. J7J 
becomes 

A T T = 0. (11) 

This equation is overdetermined and its solution is found 
by minimizing |AtT| 2 with the constraint given by Eq. 
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((5]). That is, estimation of T is given by the eigenvector 
of the smallest eigenvalue of matrix A 1 t At with a normal- 
ization factor used to satisfy Eq. ©. 

Pressure inference under an assumption of uniform tension 
(SP) 

The second method SP estimates only cell pressures. Un- 
der the assumption that all tensions are uniform, i.e., 
Tij = 1, which is strictly exact in the case of foam, Eq. 
([7]) becomes 

A P P=-A T T , (12) 

where all the components of To are 1. This problem is 
also overconditioned, and estimates of pressures are found 
by minimizing \ApP + AtTo\ 2 with respect to P. The 
solution is given as 

P = -A^AtTq, (13) 

where Ap 1 is the Moore-Penrose pseudo-inverse matrix of 
Ap. Equation (IT3l can be shown to satisfy Eq. ©. 

Bayesian inference of tensions and pressures with a prior of 
positive tension (STP) 

The third method STP is a Bayesian inference of tensions 
and pressures developed by two of us [28 . Briefly, force 
inference is carried out by MAP estimation, i.e., by taking 
the maximum value of an a posteriori distribution given 

by 

P(X) cx e -^ 2 ' 2E2 x tt(X), (14) 

where tt(X) is a prior function. We adopted the prior that 
is distributed around a positive value, because laser 
severing experiments indicate that tensions along cell con- 
tact surfaces are usually constricting in epithelial tissue 
[IIJ[24]. Hence, 

tt(X) = (2 7 :u 2 )- E ^e-' E ^ (T ^~ To)2/2uj2 x S^P,) . (15) 

A Gaussian distribution of TJy around Tq > represents 
our expectation explained above, and the Dirac S function 
is introduced to satisfy Eq. ©. We can select T = 1 by 
adjusting the scale factor. A criterion to determine hyper- 
parameters S 2 and oj 2 is to minimize the Kullback-Leibler 
distance between the probability function parameterized 
by (E 2 ,ui 2 ) (i.e., statistical model) and empirical distri- 
bution. This is formulated by maximizing marginal likeli- 
hood, or equivalently, by minimizing the Akaike Bayesian 
information criterion (ABIC) |55U5o] . 

L(S 2 ,uj 2 ) = J P(X\S 2 ,Lo 2 )dX. (16) 

We compared estimation using different priors in the arti- 
ficial data generated by a cell vertex model and confirmed 
that estimation using the prior expecting positive tensions 
gave the best fit with true values of forces. Details are de- 
scribed in [2"g] . 



2.2.3 Variable-reduction and Bayesian approaches 

It is instructive to mention how variable-reduction and 
Bayesian approaches are related. Both approaches can be 
formulated as the minimization of the function 

S(X) = \AX\ 2 +fiH(X) (17) 

with respect to X = (T, P) with constraints given by Eqs. 
(IS]) and ©. H(X) in the second term indicates our ex- 
pectation of the system, which compensates for indefinite- 
ness in the force-balance equation represented by the first 
term. Both SP (variable reduction) and STP (Bayesian) 
are formulated by using H(X) = J2[ij](Tij — T ) 2 , and 
their difference is the weight fi. 

In the variable-reduction approach, one assumes that 
the expectation given by H(X) is a strict constraint to 
be satisfied, and thus the coefficient fi is given as the La- 
grange multiplier. In contrast, in the Bayesian approach, 
S(X) is related to the posterior probability equation (Eq. 
IP)) as P(X;fx,S 2 ) cx exp (S(X; n)/2S 2 ) with /z = 
S 2 /uj 2 . The minimization of S(X) is equivalent to MAP 
estimation. The expectation H(X) is incorporated into 
the prior function, and it is not required to be strictly 
satisfied. The degree of deviation is controlled by fi, the 
weight of the second term, which is objectively deter- 
mined by maximizing the marginal likelihood L(E 2 ,ui 2 ) 
(Eq. flTBl ). By using this procedure, the first term (fit of 
data) and the second term (expectation) are balanced by 
taking into account the data quality |28) . 



3 Results 

3.1 Numerical data for foam 

Foam and an epithelial tissue are disordered cellular mate- 
rials 15,29, 38 and they have distinct mechanical natures. 
Most significantly, the tension of each contact surface is 
always kept uniform in foam |29) . Here, the geometry of a 
foam (positions of vertices and their connectivity) is ob- 
tained by numerical simulation (see Sect. 12.1.11) . and it 
is provided as an input for force inference. Although nu- 
merical simulation is carried out by considering the cur- 
vature of the membrane, we approximate it as a straight 
line when estimating forces and stress. The curvature of 
the artificial data for foam used here is not significant: 
the angular error caused by neglecting the curvature is 
5.0°(±4.9°). 

Figure [Ha) shows true bubble pressures indicated by 
a color scale. The true tensions of all contact surfaces are 

— 1. We conducted the force- inference by ST, SP, and 
STP, where the "S" stands for "straight" edges (curva- 
tures are neglected and cells are treated as polygons); "T" 
and "P" mean that tensions and pressures are unknown, 
respectively (see Sect. 12. 2^2")) . The results of force inference 
by ST, SP, and STP are shown in Figs. EKb)-(f) and are 
summarized in Table [TJ Since the pressures of all bubbles 
are set to be uniform in ST, only estimated tensions are 
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(a) true pressure (e) Estimation error of tensions 




Fig. 2. True and estimated forces for artificial foam data, (a) A color map of true pressures of the artificial data. True tensions 
are set to be T%j = 1 for all contact surfaces in the foam (not shown), (b) Tensions estimated using ST. (c) Pressures estimated 
using SP. (d) Tensions (left) and pressures (right) estimated using STP. (e) Estimation errors of tensions for ST (red) and STP 
(blue), (f) Estimated pressures plotted against true ones (top) and their errors (bottom). Pressures estimated using SP and 
STP are indicated with green and blue points, respectively. 



shown in Fig. [2Kb) . We notice that contact surfaces be- 
longing to small bubbles were estimated to have smaller 
tension. The estimation errors of tensions for individual 
contact surfaces are shown in Fig.^e) (red points). Their 
deviation from the true values (the mean residue of error) 
was a 2 = 4.5 x 1CT 2 . 

The pressures estimated using SP are indicated with 
a color map in Fig. [5Jc) , and they are plotted against the 
true ones in Fig. Hff) (green points). The results indicate 
that the accuracy of force inference was higher in SP than 
in ST (with the correlation between true and estimated 
pressures being r = 0.996 in Fig. &f), upper panel, and 
cr 2 = 3.2 x 1(T 7 in Fig. Eff), bottom panel). 

The estimated tensions and pressure obtained using 
STP are shown in Fig. [2Jd), and these estimated values 
are compared with the true values in Figs. HJe) and (f) 
(blue points). The estimation errors of tensions in STP are 
considerably smaller than those in ST (a 2 = 3.4 x 1CP 4 ). 
In addition, estimated pressures are well correlated with 



the true ones as in SP (r = 0.996; Fig. rjjf)). The error 
in the estimated pressures is very small (a 2 = 3.4 x 10~ 7 ; 
Fig. HJf), bottom panel) but slightly larger than for SP, 
where true tensions are set. 



3.2 Numerical data for the cell vertex model 

To evaluate the force- inference methods, we generate arti- 
ficial data by simulating a cell vertex model with random 
parameters (see Sect. I2.1.2[) . The results of this test for 
STP were reported previously [25] . 

Figure shows maps of true forces, and Figs. [3Jb) 
and (c) show maps of forces inferred using ST and STP, 
respectively. The pressure map obtained using SP is very 
similar to that obtained using STP (not shown). In ST, 
the estimated and true tensions do not correlate well with 
each other (r = 0.12), and the deviation of the estimated 
tensions from the true ones is a 2 = 0.17 (red points in 
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Fig. 3. True and estimated forces in artificial data obtained by cell vertex model, (a) The true tensions (left) and pressures 
(right) in the artificial data, (b) Tensions estimated using ST. (c) Tensions (left) and pressures (right) estimated using STP. 

(d) Top: Estimated tensions plotted against true ones for ST (red) and STP (blue). Bottom: Estimation errors of tensions. 

(e) Estimated pressures plotted against true ones (top) and their errors (bottom). Pressures estimated using SP and STP are 
indicated with green and blue points, respectively. 



Fig. GDJd)). SP provides good estimates of pressures with 
high correlation r — 0.97 and small mean residue of error 
cr 2 = 4.7 x 1(T 3 (green points in Fig.^e)). For STP, the 
correlations with true values of tensions and pressures are 
0.85 and 0.99, respectively (Fig. [3](d) and (e), blue points). 



the standard deviation of tensions given by ST is 0.42, 
whereas that given by STP is 0.15 (with similar results 
being obtained using all samples of the wings examined 
(n = 20)). 



3.3 Drosophila wing 

We apply the three force-inference methods to an image 
(Fig. Ufa)) of a Drosophila pupal wing (Fig. [BJc)). The 
results of estimations obtained using STP, which agree 
with laser ablation of a contact surface and the myosin 
distribution, were reported in |28) . 

Figure[3]shows the results offeree estimations obtained 
using ST (b) and STP (c). Estimated pressures obtained 
by SP exhibit maps similar to those obtained by STP (not 
shown). Since the true values of tensions or pressures are 
not known for experimental data, we compared tensions 
obtained using ST and STP and pressures obtained by 
SP and STP, respectively. The pressures obtained using 
SP and STP show a good correlation as in the artificial 
data for foam and the cell (Figs.|3Jc) and (e)); thus pres- 
sure maps obtained by the two methods are very similar 
(the pressure map obtained using SP is not shown). How- 
ever, the tensions estimated using ST show a larger de- 
viation than those estimated using STP (Figs. H[b)-(d)); 



3.4 Drosophila notum 

We applied the force-inference methods to the Drosophila 
notum (Fig. QJc) an d Fig. [S£a)). The experimental data 
used were collected at three developmental stages referred 
to as young (around 18 h APF), middle (around 22 h 
APF), and old (around 26 hr APF) (see details in [50]). 

The results of our analysis on the notum are consis- 
tent with those on the other three systems. Pressures ob- 
tained using SP and STP show a good correlation (Fig. 
GD^e)), whereas estimated tensions obtained using ST are 
more disperse than those from STP (Figs. EJb)-(d)); the 
standard deviations of tension distributions are 0.51 and 
0.11 for ST and STP, respectively (with similar results 
being obtained from all samples of the notum examined 
(n = 23)). 

Cross Validation 

The estimated stress can be cross- validated by comparing 
with that evaluated by global tissue ablation [3D], which 
indicate that the stress along the medio-lateral axis of the 
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Fig. 4. Estimated tensions and pressures for a Drosophila pupal wing, (a) An image of a Drosophila wing at 23 h APF. 
Dacatenin-TagRFP is used to highlight cell shape. Scale bar: 20 /xm. (b) and (c) Color maps of estimated tensions and 
pressures, respectively. A color scale is shown for each image, (b) Cell-junction tensions estimated using ST. (c) Tensions (top) 
and pressures (bottom) estimated using STP. (d) Comparison of estimated tensions obtained using STP and ST. (e) Comparison 
of estimated pressures obtained using STP and SP. 



notum increases during pupal development (Sect. I2.1.4|) . 
By applying the force-inference methods to images before 
ablation (the initial stage) and after the relaxation of the 
inner domain of cells (the final stage), stress tensors at 
the initial and final stages were calculated using Eq. (fTUl) 
for each sample. Since the scale factor is not determined 
by force inference, we need to reasonably calibrate it be- 
tween initial and final stages. For this, we hypothesized 
that the relationship between the cell pressure Pi and the 
cell area Ai is maintained between initial and final stages. 
In Fig. [Ha), estimated pressures are plotted against the 
cell area for initial (red) and final (blue) stages. The fit- 
ting function P(A) = a/^/A + b is given by a dimensional 
argument considering Laplace's law, as in the case of foam 
[29] . By evaluating coefficients a and b by the least-squares 
method, we obtained the fitting curves indicated by dot- 
ted lines in Fig.jBJa). Then, the scale is calibrated for these 
two lines (Fig. [6jb)). 

Stress tensors for initial and final stages, cr l and cr* , 
are estimated using STP at distinct developmental stages. 
We measured the normal stress difference a a = {p yy — 
o~xx)/2 and a xy in Fig. [7]Ja), and calculated the differ- 
ence of these quantities between the initial and final stages 
[Acta = (J a~ (J a an( l Ao X y = a xy — al y ), because they 
are independent of the unknown additive constant in the 
pressure, and characterize change of mechanical state in- 
duced by the ablation. The estimated amplitude and the 
difference of o~ xy were smaller than those of a a in older pu- 
pae, indicating that the stress in the notum was stronger 
along the y axis (medio-lateral axis) (Fig. E£a)). However, 



a A at young and middle stages exhibit weaker changes 
upon ablation of cells. These results qualitatively agree 
with those obtained from global ablation [30]. To quanti- 
tatively cross-validate the two methods, Ao a values are 
directly compared (Fig. E£b)). Aoa values obtained by 
global ablation and STP exhibit a good correlation (with 
a correlation coefficient of r = 0.64 for calibrated data 
and r = 0.59 for uncalibrated data). By repeating the 
same procedures in SP and ST, we obtained r = 0.63 and 
r = 0.59 for calibrated and uncalibrated data using SP, 
and r = 0.56 for uncalibrated data using ST. 

3.5 Robustness to image processing error 

We evaluated robustness of the force-inference methods 
to image processing error. In [25], a similar test for STP 
using wing data was briefly mentioned. We made 100 sam- 
ples by adding Gaussian noise to position coordinates of 
vertices in the original data, where the standard devia- 
tion of the noise is 5% of the mean length of cell con- 
tact surfaces. Then we checked how estimated tensions 
and pressures deviate among the noised data. Tests were 
conducted for data from the foam simulation, vertex cell 
model, and Drosophila notum. For reference, among the 
sequential images of notum taken in a very short period of 
time (about 1.6 seconds), the mean deviation of vertices 
positions extracted by image processing was less than 3% 
of the mean contact surface lengths. 

The deviations of estimated tensions among samples 
are large for ST; mean relative values of standard devi- 
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Fig. 6. The calibration of scale factors of the stress tensor, 
(a) Estimated values of cell pressures plotted against cell area 
before cells were ablated (red) and after the tissue is relaxed 
(blue), (b) The area-pressure relation is assumed to be main- 
tained between the two time points. By using a fitting function 
P{A) = a/VA+b, the pressures at the initial stage can be cal- 
ibrated to coincide with those at the final stage by selecting a 
scale factor a and by adding a hydrostatic value of pressure b. 
Only the scale factor is used for the following analysis. 




-0.2 -0.1 0.1 0.2 




(d). 



(e) 



a 

to 
E 

HI 




Q_ 




Estimated tenstion by STP 



-0.1 0.1 0.2 

Estimated pressure by STP 



Fig. 5. Estimated tensions and pressures for a Drosophila 
pupal notum. Data from [3DJ are used, (a) An image of a 
Drosophila notum. Scale bar: 10 /im. (b) and (c) Color maps 
of estimated tensions and pressures. A color scale is shown for 
each image, (b) Cell-junction tensions estimated using ST. (c) 
Tensions (top) and pressures (bottom) estimated using STP. 
(d) Comparison of estimated tensions obtained using STP and 
ST. (e) Comparison of estimated pressures obtained using STP 
and SP. 



ations (standard deviation divided by mean tension) are 
17%, 69%, and 73% for foam, cell vertex model, and no- 
tum, respectively. On the other hand, those for estimated 
tensions by STP were much smaller: 5.0% (foam), 7.8% 
(cell vertex model), and 8.5% (notum). Deviations of es- 
timated pressures for SP were 9.5 x 10~ 4 , 0.16, and 0.020 
for foam, cell vertex model, and notum, and those for STP 
were almost the same: 1.5 x 10~ 3 , 0.19, and 0.026. These 
values are sufficiently small compared to the estimated 
dispersion of individual cell pressures, indicating highly 
robust pressure estimations in SP and STP. 
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Fig. 7. Stress inference in the notum. (a) Difference of es- 
timated stresses obtained using STP between initial (before 
laser ablation) and final (after the relaxation) stages (Acta = 
CT 4 — a A an d Aa xy = al y — o\.y). (b) Estimated Aoa values 
obtained using STP plotted against those obtained by global 
ablation of tissue in [3DJ . Different colors indicate developmen- 
tal stages of samples in all panels (yellow: young, red: middle, 
and blue: old). 



Moreover, we found that maps of tensions obtained 
by ST showed "patches" (distinct regions where the ten- 
sion seems locally uniform) and that positions of these 
patches differed among individual noised samples (Fig. [8]). 
These patches were also seen in tension maps estimated 
by ST from the original data of the wing and notum (Fig. 
|4jb) and Fig. [5fb)). The large error in the estimated ten- 
sions using ST can be explained by the appearance of the 
patches that are sensitive to errors in vertices positions. 

Finally, the robustness in the stress inference was ex- 
amined. The deviation among noised samples is largest in 
ST: for foam, the deviations of cta are 1.1 x 10~ 4 in ST, 
5.8 x 10~ 5 in SP, 8.8 x 10~ 5 in STP. For the cell vertex 
model, they were 1.0 x 10~ 2 , 3.6 x 10~ 3 , and 5.7 x 10~ 3 in 
ST, SP, and STP, respectively. The difference between ST 
and STP was more significant in the experimental data 
of the wing (Fig. 01a)), where the deviations of a a were 
1.3 x 10~ 2 in ST, 5.1 x 10" 4 in SP, and 1.1 x 10~ 3 in STP. 
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Fig. 8. Examples of estimated tensions by ST for data with 
errors in the vertex positions. Two samples are shown for (a) 
simulated foam data, (b) a simulated cell population from the 
cell vertex model, and (c) Drosophila notum (the method to 
obtain them is described in Sect. l3.5j) . The estimated values of 
tensions and the position of patches vary among samples. 

For the sample of notum shown in Fig. [5ja), they were 
5.5 x 1(T 3 in ST, 4.9 x 1CT 4 in SP, and 8.4 x 1CT 4 in STP. 



4 Discussion 

4.1 Force-inference methods 

In the present study, we performed a comparative anal- 
ysis of force and stress inference in tissue. We employed 
three types of force-inference methods, which require dif- 
ferent assumptions on the unknown variables (tension and 
pressure) for treating indcfinitcncss in the inverse problem 
between forces and cell shape; ST (resp.: SP) estimates 
only tensions (resp.: pressures) under the assumption that 
pressures (resp.: tensions) are uniform, and STP estimates 
both by having a prior that tensions are positive. We pre- 
pared four different data sets, in which the assumptions 
in each force-inference method are either: strictly exact, 
reasonable, incorrect, or not checked, which enables us to 
better evaluate the force-inference methods (see also Table 

The assumption in SP that all tensions are uniform is 
strictly exact in foam and is not exact in the cell vertex 
model or in tissue (and the variance of tensions may differ 
in each tissue). Estimates of pressures exhibited high ac- 
curacy in both the numerical data from foam and the cell 
vertex model (Table [1]), suggesting that the difference in 
the cell area is a good indication of cell pressures in the ar- 
tificial data employed in this study. Whether SP provides 
reasonable estimates of pressures in tissues where tensions 
vary greatly among contact surfaces is a subject for future 
study. 

The assumption in ST that all pressures are uniform 
is incorrect in foam and in the cell vertex model and is 
not directly checked in tissue. The errors of tension infer- 
ence obtained using ST are large in both the numerical 
data from foam and the cell vertex model. In foam, the 
estimated tensions exhibit a positive correlation with the 



length of the contact surface, which disagrees with the 
physical nature of foam. In the cell vertex model, esti- 
mated tensions are significantly more disperse than the 
true ones. 

The large estimation error in ST may result from its in- 
correct assumption. Another, but not mutually exclusive, 
possibility is indicated by our results: adding noise in in- 
put data causes the appearance of " patches" in a map of 
tensions. Note that the force-balance equations © and @ 
provide only local information and that the force-balance 
equations become non-exact upon addition of error in the 
vertex positions. Thus, one can speculate that a gradual 
spatial change in tensions and pressures (i.e., modes with 
long wavelength) can be allowed, which generates a spa- 
tial modulation comparable to system-size. The patches 
of tensions generated by long-wavelength modes are sen- 
sitive to image processing error as indicated by our ob- 
servation that the positions of the patches vary among 
noised samples (Fig. [BJ). The positions of the vertices and 
the force-balance equations are exact in the original arti- 
ficial data, while image processing errors are unavoidable 
in experimental data. This is the reason why the patches 
of tensions were much more evident in the experimental 
data than in the artificial ones (compare Fig. [2jb) and 
Fig. Mb) with Fig. \Mb) and Fig. EJb)). In contrast, the 
force inference by STP is more accurate and more robust 
to image processing error. In fact, in STP, the prior that 
expects Tij should be close to To > works as the 'reg- 
ularization term' to avoid over-fitting by eliminating the 
long-wavelength mode, which makes the force inference 
more robust (see also Sect. 12.2.31) . 

The assumption in STP that tensions are distributed 
around a positive value is reasonable in epithelial tissue, 
as suggested by laser ablation of individual contact sur- 
faces [TTTI24] . Our results clearly indicate that estimates 
of pressures and tensions correlate well to the true values 
in both the numerical data from foam and the cell vertex 
model. To summarize, STP gives better estimates of ten- 
sions than ST. SP is slightly better at pressure inference 
than STP for foam data (Fig. H£f)), in which the assump- 
tion = 1 is the true physical nature of the system, 
whereas STP gives slightly better results than SP for the 
data generated by the cell vertex model, in which tensions 
are not uniform (Fig. [3Jf)). 

The application of force-inference methods to exper- 
imental data yielded results consistent with those from 
artificial data. Estimates of pressures by SP and STP are 
highly correlated in the Drosophila wing and notum. In 
contrast, estimates of tensions by ST vary more widely 
than those by STP; several contact surfaces are estimated 
to have almost zero tension, which may result from the 
spatial distribution of tensions where the strength of ten- 
sions are significantly different among neighboring patches. 

4.2 Summary of the cross-validation 

We compared two different stress measurement methods 
in tissue: One is robust and model-independent [3D] and 
one is non-invasive and can yield a space-time map (i.e., 
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Table 1. Estimation errors of force-inference methods for artificial data. The errors are represented by the mean residues a 2 . 



Foam simulation data Cell vertex simulation data 
estimation error ST SP STP ST SP STP 

tension 4.5 x 1(T 2 3.4 x 10" 4 1.7 x 1CT 1 3.8 x 10" 3 

pressure 3.2 x 10" 7 3.4 x 10" 7 4.7 x 10" 3 1.0 x 10~ 3 



STP) [55] • Annular ablation of tissue showed that the 
stress of the notum increases along the medio-lateral axis 
during pupal development. STP can detect such devel- 
opmental changes in stress. Moreover, the value of Aua 
inferred using STP better agrees with experimental data 
than those found using ST and SP. These results serve as 
double-checks for the anisotropy of stress in the notum 
and in turn this reinforces the validation of the annular 
ablation method of stress measurement. 



4.3 The curvature of a contact surface 

In the present study, we neglected the curvature of the cell 
contact surface and obtained good estimates of pressures 
in SP and STP. Note that neglecting the curvature does 
not imply that the difference in pressures among cells is ig- 
nored. It was previously discussed that the error of stress 
evaluation under the straight-edge approximation could 
be small, if the curvature is small [35] . Neglecting the cur- 
vature was also demonstrated to give a good estimation in 
vivo by checking Laplace's law for the tension estimated 
under the straight-edge approximation [27] . One advan- 
tage of neglecting the curvature of the contact surface is 
that we can bypass having to make a precise curvature 
measurement from an image, which is known to be dif- 
ficult and/or error-prone |40j . If the pressure difference 
between cells is large, the curvature cannot be neglected; 
hence the force-inference method needs to be extended to 
use information about the relationship between forces and 
the curvature. 



4.4 A practical choice of force and stress inference 

Selection of the most suitable force-inference method is 
dependent on the nature of the system and experimen- 
tal design. As summarized above (Table [1} , tension in- 
ference obtained using STP is better than that obtained 
using ST in accuracy and in robustness and STP gives 
more robust stress inference than ST, whereas pressure 
and stress inference obtained using STP and SP are com- 
parable. On the other hand, SP and ST involve solving 
a linear equation only once, STP requires performing QR 
decomposition many times during the maximization pro- 
cess 28 . Although it requires less than ~10 min for most 
of the data investigated in this study (~200 cells), com- 
putational time increases for data with a larger number of 
cells. Collectively, STP should be the first choice for ten- 
sion, pressure, and stress inference. SP can be useful for 
stress and pressure inference, when the variance of ten- 
sions are known to be sufficiently small as in the notum, 



and when there are large numbers of cells (e.g., more than 
thousands) in the system of interest. 

In conclusion, the present study strengthens the va- 
lidity of our force-inference [55] and stress measurement 
30 methods. Their future improvement would accelerate 
studies of the physical regulations of development. 
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